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Abstract 



We use the invariant mass distribution of Drell-Yan dileptons as measured by 
the CDF and D0 Collaborations at the Fermilab Tevatron and make a careful 
analysis to constrain Kaluza-Klein models with large extra dimensions. The 
combined data from both collaborations lead to a conservative lower bound on 
the string scale Ms of about 1 TeV at 95% confidence level. 
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Recently, the idea that gravity could become strong at scales of the order of 
a few TeV has attracted a great deal of attention |T|]. This is made possible if 
we allow for large compactified dimensions at the TeV scale. While such ideas 
can be fitted in within the scheme of quantum field theories]^], a more natural 
construction p, |J involves string theories with all Standard Model (SM) fields 
living on a three-dimensional D-brane (or 3-brane) embedded in a space of (4 + d) 
dimensions (bulk). Of course, the original suggestion that we live in a spacetime 
continuum with more than the three canonical spatial dimensions was made early 
in this century ||, but these Kaluza-Klein (KK) theories, as they are called, have 
not been able to satisfactorily reproduce the observed mass spectrum. Such 
ideas, however, have always formed a basic ingredient of string theories^]. In 
fact, models having extra dimensions with compactification scales of the order of 
a few TeV have been proposed |7j from time to time in the literature with various 
motivations. However, it is the discovery of D-branes|| which has provided the 
rather venerable KK theories with a new lease of life over the past year. 

In a nutshell, the ideas proposed by Arkani-Hamed, Dimopoulos and Dvali 
(ADD)[[| and by Antoniadis et a/. [[J follows. They suggest — as all KK 

theories do — that spacetime consists of (4 + d) dimensions. The extra (spatial) 
d dimensions are compactified, typically on a (^-dimensional torus T d with radius 
R each way. Since gravity experiments have not really probed the sub-millimetre 
regime, it is proposed that R can be as large as ~ 0.1 — 1 mm, a very large value 
when compared with the Planck length ~ 10~ 33 cm. Though the actual value 
of Newton's constant G% in the bulk is of the same order as the electroweak 
coupling, its value G^ in the effective 4-dimensional space at length scales ^> R 
is the extremely small one measured in gravity experiments. This is decribed by 
a simple relation derived 0] from Gauss' Law, 



4? 



2 r -i (d+2) 

~ R d 



where Mpi ~ 1/ y/G^ denotes the Planck mass. If Mp^^ ~ 1 TeV, then R ~ 
10 3o/d-i9 m . This means that for d — 1, R ~ 10 11 m, which, in turn, means that 
deviations from Einstein gravity would occur at solar system scales; since these 
have not been seen, we are constrained to take d > 2. For these values R < 1 mm, 
hence there is no conflict with known facts. It is also perhaps worth mentioning 
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that we would normally require d < 7, since that is the largest number allowed 
if the string theory is derivable from M-theory, believed to be the fundamental 
theory of all interactions. In the ADD model the smallness of Newton's constant 
is a direct consequence of the compactification-with-large-radius hypothesis and 
hence there is no hierarchy problem in this theory^. 

In traditional KK theories, the mass-spectrum of nonzero KK modes arising 
from compactification of fields living in the bulk is driven to the Planck scale 
Mpi . This problem is avoided in the ADD model by having the SM particles live 
on a 'surface' with negligible width in the extra d dimensions, which we identify 
with the 3-brane. The SM particles may then be thought of as excitations of open 
strings whose ends terminate on the brane; gravitons correspond to excitations of 
closed strings propagating in the bulk. Thus, the only interactions which go out 
of the 3-brane into the bulk are gravitational ones. We thus have a picture of a 
4-dimensional 'surface' embedded in a (4 + d) -dimensional space, where SM fields 
live on the 'surface', but gravitons can be radiated-off into the bulk. Noting that 
the SM fields are confined to the 3-brane, it is obvious that the only new effects 
will be those due to exchange of gravitons between particles on the 3-brane. To 
construct an effective theory in 4 dimensions, gravity is quantized in the usual 
way, taking the weak- field limit, assuming that the underlying string theory takes 
care of ultraviolet problems. The interactions of gravitons now follow from the 
(4 + d) -dimensional Einstein equations in the compactification limit. Feynman 
rules for this effective theory have been worked out in detail in Refs. [Ill] and 



[12]. We use their prescriptions in our work. On the 3-brane, the couplings of 
the gravitons to the SM particles will be suppressed, as is well-known, by the 
Planck scale Mpi ~ 1.2 x 10 19 GeV. This is offset, however, by the fact that, 
after compactification, the density of massive KK graviton states in the effective 
theory is very high, being, indeed proportional to Mpi /M^ + . The Planck 
mass dependence cancels out, therefore, leaving a suppression by the string scale 
Ms = M p A ^ d ^ ~ Mew- in the ADD theory, therefore, the tower of KK graviton 
states leads to effective interactions of electroweak strength. A further assumption 
made in our work — and in other phenomenological studies — is that F-particles, 



5 A related problem, that of stabilization of the compactification scale, exists, however; this 
has been discussed in Ref. 01. 
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excitation modes of the 3-brane itself in the bulk — are heavy and do not affect 
the processes under consideration. This corresponds to a static approximation 
for the brane. It is also relevant to mention that the dilaton field associated with 
the graviton couples only to the trace of the energy-momentum tensor, i.e. to 
the mass of the SM particles at the vertex. For light fermions, as we have in the 
Drell-Yan process, this means that the interactions of the dilaton can be safely 
neglected. 

Using these Feynman rules, it has been possible to explore a number of dif- 
ferent processes where the new interactions could cause observable deviations 
from the SM. Only two new parameters enter the theory: one is the string scale 
M s = M ( p+ d) . The other is a factor A, of order unity and indeterminate sign, 
which arises when we sum over all possible KK modes of the graviton. As the 
amplitudes for virtual graviton exchange (with which we are concerned in this 
work) are always proportional to A/M|, it is usual to absorb the magnitude of 
A into Ms; this reduces the uncertainty to A = ±1. Obviously this determines 
whether the graviton exchanges interfere constructively or destructively with the 
SM interactions. 

Remembering that the gravitons couple to any particle with a non-vanishing 
energy-momentum tensor, it is possible to make a variety of phenomenoogical 
studies of the new interactions and to test the workability of the ADD model. 
Though the phenomenology of this model has not yet been fully explored, several 
important results are already available in the literature. These can be classi- 
fied into two types: those involving real KK graviton production, and those 
involving virtual graviton exchange. A real KK mode of the graviton will have 
interactions with matter suppressed by the Planck scale M$ and will therefore 
escape the detector. One can, therefore, see signals with large missing momen- 
tum and energy if an observable particle is produced in association with a KK 
graviton mode. However, cross-sections for these depend explicitly on d, the num- 
ber of extra dimensions, and bounds derived from data reflect this dependence. 
Some of the processes examined so far include single-photon final states at e + e~ 



colliders [11, 13, 14 1 as well as hadron colliders [11], monojet production at hadron 
colliders fTT|, 13]] , two-photon processes at e + e~ colliders^], single-Z production 
at e + e~ colliders |14[] and the neutrino flux from the supernova SN1987A[|1(], |I~6|[ . 
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Each process can be used to obtain a bound on the string scale Ms for a given 
number d. The most dramatic of these bounds is Ms > 50 TeV for d = 2 and it 
comes from a study|l| of neutrinos from the supernova SN1987A. However, this 
last bound drops to about a TeV as soon as we go to d = 3. Most of the other 
processes lead to lower bounds of about 1-1.1 TeV on the string scale for d = 2, 
but these bounds become much weaker for d > 3. 

Virtual (KK) graviton exchanges lead to extra contributions to processes in- 
volving SM particles in the final state and can be observed as deviations in the 
cross-sections and distributions of these from the SM prediction. After summa- 
tion over all the KK modes of the graviton, the final result is proportional to 
sgn(A)/Mj, with practically no dependence on the number of extra dimension^ 
Each process can be used to obtain a bound on the string scale Ms for a given 
sign of A. Some of the processes examined include Bhabha and M /Slier scatter- 



ing at e + e~ colliders pi], Qn photon pair-production in e + e~[TT], [22| and hadron 
colliders [TT], fermion pair production in 77 colliders ||17||, Drell-Yan production 
of dileptons[T8H, dijet[19] and top-quark p0| pair production at hadron colliders, 
deep inelastic scattering at HERA^T], [17]], massive vector-boson pair produc- 
tion in e + e~ collisions |22|, |23j and pair production of scalars (Higgs bosons and 
squarks) at both e + e~ and 77 colliders [[H]. Among the best of these bounds is 
M s > 920 (980) GeV for A = +1(-1) which comes from a study|fr8[ of exper- 
imental data on Drell-Yan leptons at the Tevatron. We make a more elaborate 
analysis if the same data in this work. 

The contributions to the Drell-Yan production of dileptons at hadron collid- 
ers from graviton exchanges have been considered by Hewett JTBj] . Some of her 
findings relevant to the Tevatron are: 

• There is very little difference between the cases A = +1 and A = — 1 for the 
dilepton invariant mass distribution. 

• The A = ±1 cases differ, however, in the angular distribution; therefore, 
widely differing forward-backward asymmetries may be predicted. 



6 This is really because the density of graviton KK modes is approximated by a continuum, 
as a result of which mass degeneracies due to the number of extra dimensions are lost, at least 
to the leading order. In a sense, therefore, bounds from virtual graviton exchange are more 
general. 
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• There are large deviations between the SM and the ADD model for large 
invariant masses. 

• The gluon-gluon contribution to the Drell-Yan process (see below) is much 
suppressed compared to the quark-initiated process. 

• The bounds can increase to about 1.15 (1.35) TeV for A = +1(— 1) in Run-II 
of the Tevatron. 

We agree with most of these results at the generator level. However, in the 
absence of published details about the angular distribution of dileptons observed 
by the CDF and D0 Collaborations, we confine our analysis to the invariant mass 
distributions only. Hence we do not make a separate analysis for the two signs 
of A. 




Figure 1. Feynman diagrams for the contribution to the Drell-Yan process from (a) the 
Standard Model and (b,c) exchange of a Kaluza-Klein graviton. 

The Drell-Yan cross-section, including the effects of Kaluza-Klein graviton 
exchanges, is given by the above Feynman diagrams. The Standard Model dia- 
grams (a) involving exchange of a photon or a Z-boson in the s-channel, interfere 
with the diagram with s-channel exchange of a Kaluza-Klein graviton (b), while 
the diagram (c) has no Standard Model analogue. 

Evaluating these leads to the result 

ct D y(pP -> ^t~) = J dxidx 2 fg/pfa) f g /p(x 2 ) o(gg -> (1) 
+ Yl J dxidx 2 [f q / P (xi) fq/p(x 2 ) + fq/p(xi) f q /p{x 2 )\ v(qq -> £ + T) , 

q=u,d,s 
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where f a /b(x) denotes the flux of a parton a in a beam of particles b, 

1 



a(qq or gg -> tt ) 



167TS 2 



\M(qq or — > 



-M2 



(2) 



and \M.\ 2 represents the squared Feynman amplitude summed over final spins 
and averaged over initial spins and colours. 

Evaluation of the Feynman diagrams gives, for the gluon-induced process 
(which has no Standard Model analogue): 

2 r 

s 4 + 2tu (t - uf 



(3) 



when all the graviton Kaluza-Klein modes have been summed over. 

Evaluation of the Feynman diagrams gives for the quark-induced process (in- 
cluding interference terms): 



\M( q q^en-)\z = Tii + T?? K , 



(4) 



where the Standard Model contribution is given below. We adopt the convention 
that T a denotes the contribution from exchange of a particle a and T ab denotes 
the interference term between diagrams with exchange of a and b respectively. 
With these, we get 
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L e = 4 sin 2 9 W - 2 , = 4 sin 2 W , 

= 4(T 3<Z - Q q sin 2 W ) , i? g = -4Q q sin 2 ^ , 
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for the couplings. The non-Standard part, using the same convention, is given 

by 
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when, as before, all the Kaluza-Klein modes have been summed over. 

The above formulae represent the lowest order (LO) calculation in perturba- 
tion theory. The calculation of higher-order effects, especially next-to-leading or- 
der (NLO) and next-to-NLO (NNLO) QCD corrections has been done in detail p6 



for the SM process represented by Tg q M . No corresponding calculations have been 
attempted as yet for the KK parts, T% K and M.{gg — > In the absence 

of such a calculation, we make the assumption that the change in the LO cross- 
section due to QCD corrections — the 'K-factor' — is identical for the SM and 
KK parts. Our results are, therefore, correct only within this approximation^. 
However, we do not expect a proper calculation of NLO effects to make a dras- 
tic change in our rough-and-ready results, because the dominant contribution to 
dilepton production at the Tevatron comes from quark-induced processes. Since 
the SM and KK results both arise from colour-singlet exchange, the actual 'K- 
factor' is likely to be rather similar in both cases. For gluons, this is not true, but 
the gluon-induced process makes only a minor contribution at Tevatron energies. 

In keeping with this philosophy, therefore, we have extracted, for each value 
of the dilepton invariant mass M = M1+1-, a 'K-factor' by taking the ratio of the 
LO SM cross-section calculated using the above formulae with that calculated 



using the full NNLO calculation of Ref. f26j. This set of ratios is then used to 
scale the entire differential cross-section when the KK effects are included. It 



7 This places our results on an equal footing with a large number of experimental bounds on 
new physics scenarios, such as those involving quark and lepton compositeness[^, for which 
the QCD corrections are not available. 
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is worth pointing out that this procedure also takes care of the leading effects 
arising from initial-state radiation. Finally, it is relevant to mention that we have 



used the CTEQ-4M set of structure functions p5[| to calculate the initial state 
parton luminosities. 

We now describe our analysis in some detail. The D0 Collaboration has 
presented [f27| the e + e~ invariant mass distribution in 9 bins starting from 120 



GeV till 1 TeV using the di-electron data collected with 120 pb _1 of luminosity. 
The cuts relevant for the cross-section calculation are given below. No distinction 
is made between the electron and the positron. 

• The transverse momentum of both the isolated electrons must satisfy pr > 
25 GeV. 

• The electrons are called CC (for Central Calorimeter) if they satisfy \r]\ < 
1.1, 77 being the pseudorapidity; they are called EC (for End Cap) if they 
satisfy 1.5 < \r)\ < 2.5. 

Only those events are considered in which there is at least one CC electron, while 
the other can be CC or EC. The acceptances described above are taken into 
account while estimating our Monte Carlo cross-sections. These cross-sections 
need to be further convoluted with efficiencies []2?| which are (74.1 ± 0.6)% when 



both electrons are CC and (52.6 ± 1.0)% when one of them is EC. Multiplying 
by the luminosity now gives us a prediction for the number of di-electron events 
expected in each mass bin, which is then compared with the D0 data. 

The CDF Collaboration has presented |28| results for dimuon samples, using 



107 pb of data. The relevant cuts are given below. 

• The reconstructed rapidity y of the virtual s-channel state ('boson rapidity') 
is required to satisfy \y\ < 1 for all events. 

• Both muons are required to satisfy \r]\ < 1, which confines the analysis to 
the central region. 

• A back-to-back cut \r]i + 1] 2 \ > 0.2 is imposed: this gets rid of cosmic ray 
backgrounds. 



S 



• Both muons are required to satisfy a 'loose' transverse momentum cut of 
Pt > 17 GeV and at least one is required to satisfy a 'tight' cut of pr > 20 
GeV. 

These cuts are applied in our Monte Carlo generator to estimate the cross-section 
times acceptance for the 6 mass bins in the range 120 GeV to 500 GeV presented 
in Ref. [[^J (Table X). These are convoluted with the experimental efficiencies 
(Table VI of Ref. f28fl). We then obtain an additional correction factor for each 



mass bin by normalising our SM expectation to the numbers given in Ref. [28 



This may be expected to take care of the effect of other detector-specific cuts like 
triggers, etc. Finally, we use this correction factor along with our generator-level 
acceptance and the experimental efficiencies to estimate the number of events in 
each mass bin for various values of Ms- The choice of only 6 mass bins in the 
range 120 GeV to 500 GeV is because the ADD model predicts wider deviations 
from the SM in the higher mass bins (see Fig. 2). We also take note of the fact 
that no events are seen at CDF in the mass bin 500 GeV to 1 TeV. 
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Figure 2. Illustrating the effects of TeV scale quantum gravity on the invariant mass distri- 
butions of dileptons seen at the Tevatron by the D0 and CDF Collaborations respectively. Solid 
lines show the SM prediction; dashed lines show the predictions of the ADD model for marked 
values of Ms ■ 

In Fig. 2 we show the differential cross-section as a function of the invariant 
mass M of the dilepton, compared to the D0 and CDF data. We have set A = +1, 
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but A = — 1 will not make a discernible change in the figure. Solid lines show the 
SM prediction; dashed lines show the predictions of the ADD model for M s = 0.5, 
0.75, 1, 1.25 and 1.5 TeV respectively. The data points correspond to those used 
in our analysis and do not represent the full set of available points. Error bars 
are presented at 68% confidence level (C.L.) if there are events in the relevant 
mass bin and a 95% C.L. upper bound if there are no events in the relevant mass 
bin. The D0 numbers correspond to a differential cross-section da/dM: this 
is obtained by dividing the cross-section (modulo cuts) in the mass bin by the 
width of that bin. The CDF numbers correspond to a double differential cross- 
section d 2 a j dMdy in both M and y: this is obtained by dividing the cross-section 
(modulo cuts) in the mass bin by the width of that bin as well as by a factor 



As is apparent from the figure, the string scale cannot be anywhere near 500 
GeV, since that would show extreme deviations from the observed data. This 
is just one of the arguments which tells us that quantum gravity effects must 
lie at scales of a TeV or more. On the other hand, as Ms approaches 1 TeV, 
the differentiation between signal and background is less striking. This is partly 
because the deviations arise only in the high mass bins, where no events are 
expected with the current luminosities. 

The actual limits on the string scale M s are calculated using a Bayesian 
analysis of the shape of the mass distribution of events. For a value M s of the 
string scale, the expected number of events in the k th mass bin can be written 
as: 



where C is the data luminosity, b k is the expected background, e k is the dilepton 
detection efficiency and cr k (Ms) is the expected dielectron cross section with 
inclusion of the effect due to large extra dimension. 

The posterior probability density for the string scale to be Ms, given the 
observed data distribution (D), is given by 



Ay = 2. 



N k (M s )=b k + Ce k a k (M s ) 



(8) 



P(M S \D) = 1/ 



db de dC Y\ 
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n 



e - Nk ( M s)N k (M s ) 



P(b,Ce) P(M S ). 



(9) 



10 



In the above equation the term in square brackets is the likelihood for the data 
distribution to be from a model with string scale M s . The prior probability 
P(b, He) is taken to be a product of independent Gaussian distributions in b, C 
and e, with the measured value in each bin defining the mean and the uncertainity 
defining the width. The overall factor l/A is just a normalisation. Since the 
excess cross-sections due to graviton exchanges are combinations of direct terms 
proportional to 1/Mj and interference terms proportional to 1/Mj, we consider 
a prior distribution P(Ms) uniform in (a) 1/Mj and (6) 1/M| separately. The 
limit on the string scale from a prior uniform in 1/Mj represents a conservative 
estimate; using a prior uniform in 1/Mj provides more stringent limits. From the 
above posterior probability, the cumulative probability = J^ s P(M' s \D)dM' s can 
be calculated. The Ms value at which the cumulative probability equals 0.95 is, 
then, the 95% C.L. limit. We also combine the data using the simple expedient 
of treating the CDF probability as a prior for the DO analysis (and vice versa). 
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Figure 3. Showing the (cumulative) posterior probability for the ADD model with different 
values of the string scale Ms assuming a prior probability which is (a) uniform in 1/Mj and 
(b) uniform in 1/Mj. 

In Fig. 3, we have plotted the cumulative posterior probability for the ADD 
model with string scale Ms as a function of (a) 1/M| and (b) 1/Mf. The dashed 
(dash-dot) lines indicate the results of considering the D0 (CDF) data alone, 
while the solid lines show the result of a combined fit. A glance at the figure will 
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show that the horizontal (dotted) lines correspond to 95% CL. limits, while the 
fact that the curve saturates for higher values of l/M^' 8 ^ shows that the SM values 
(Mg — > oo) constitute the best hypothesis to fit the data. In a more quantitative 
idiom, we may interpret the vertical (dotted) lines as lower bounds on the string 
scale Ms- It is clear that the bound of 927 GeV assuming a prior probability 
of l/Mg, and using the CDF data, is consistent with that reported in Ref. 



while the value 874 GeV assuming a prior probability of 1/ Mf represents a more 
conservative estimate. The DO data provide an improvement^ in the bound by 
about 100 GeV in both cases. Since the cross-section varies principally as 1/Mf 
in the region around Ms = 1 TeV (this is reflected in the fact that it depends 
very weakly on the sign of A), this corresponds to an increase in the sensitivity 
by a factor of about 2.6. Combining the data increases the sensitivity by another 
factor of about 1.2, which takes the bound to 1080 (1016) GeV, depending on the 
choice of prior probability. Increasing the energy to 2 TeV and the luminosity to 
2 fb _1 , which may be expected with the commissioning of the Main Injector in 
Run-II, improves the bounds by a further 200-300 GeV; this corresponds to an 
improvement in the sensitivity by a factor close to 4. 

To conclude then, we have used published dilepton data from the D0 and CDF 
Collaborations to put bounds on the string scale Ms- This is the fundamental 
scale of the ADD model, which envisages large compact dimensions in addition 
to the known (noncompact) ones and predicts strong quantum gravity effects 
at TeV scales. Only the invariant mass distribution has been used and not the 
angular distribution. The latter might show some sensitivity to the sign of A. For 
the current analysis, however, there is hardly any such sensitivity. Our result is 
also independent of the number of extra dimensions d. We obtain a bound on 
M s of 900 GeV (900 GeV - 1 TeV) using CDF (DO) data alone and a bound 
of around 1.0 - 1.1 TeV using the combined data from both experiments. This 
is one of the most stringent bound obtained from collider studies at the present 
time and is likely to be improved (to about 1.3 TeV) in Run-II of the Tevatron. 

The authors would like to thank Ashoke Sen and K. Sridhar for reading 

8 This is for two reasons: (a) the published results from D0 use a slightly higher integrated 
luminosity and (b) the D0 Collaboration presents more data in the higher mass bins — where 
most of the deviations lie — than the CDF Collaboration. 
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References 

[1] See T. Banks, M. Dine and A.E. Nelson, SCIPP preprint, SCIPP-99-03 (Mar 1999), |hcp~- 
th/9903019| , for a lucid discussion of the theoretical issues involved. 

[2] N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B429, 263 (1998) . 

[3] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B463, 257 
(1998) . 

[4] R. Sundrum, Phys. Rev. D59, 085009 (1999) and Phys. Rev. D59, 085010 (1999) . 

[5] G. Nordstrom, Phys. Zeitschr. 15, 504 (1914); Th. Kaluza, Situngober. Preuss. Akad. 
Wiss. Berlin, p.966 (1921); O. Klein, Z. Phys. 37, 895 (1926), 1938 Conference on 
New Theories in Physics, Kazimierz, Poland (1938) and Helv. Phys. Acta, Supp.IV, 58 
(1956); For a fascinating discussion of the early history of gauge and KK theories, see 



L. O'Raifeartaigh and N. Straumann, hcp-ph/9810524; a more technical review is given 



by D. Bailin and A. Love, Rept. Prog. Phys. 50, 1087 (1987). 

[6] M.B. Green, J. Schwarz and E. Witten, Superstring Theory, Vols. 1 and 2 (Cambridge 
University Press, 1987). 

[7] V. Rubakov and M. Shaposhnikov, Phys. Lett. B125, 136 (1984) ; A. Barnaveli and 
O. Kancheli, Sov. J. Nucl. Phys. 51, 573 (1990); I. Antoniadis, C. Muhoz and M. Quiros, 
Nucl. Phys. B397, 515 (1993) ; I. Antoniadis, K. Benakli and M. Quiros, Phys. Lett. 
B331, 313 (1994) . 

[8] For an excellent pedagogical review, see, for example, J. Polchinski, Tasi Lectures on 
D-Branes, |hcp-th/9611050| (1996). 



[9] N. Arkani-Hamed, S. Dimopoulos and J. March-Russell, SLAC preprint SLAC-PUB-7949 



(1998), |hcp-th/9809124 . 



[10] N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Rev. D59, 086004 (1999) 
[11] G.F. Giudice, R. Rattazzi and J.D. Wells, Nucl. Phys. B544, 3 (1998) . 



[12] T. Han, J.D. Lykken and R.-J. Zhang, Madison preprint MADPH-98-1092 (1998), hep 



ph/9811350 



13 



[13] E.A. Mirabelli, M. Perelstein and M.E. Peskin, SLAC preprint SLAC-PUB-8002 (1998), 
tiep-ph/9811337 . 



[14] K. Cheung and W.-Y. Keung, Davis preprint UCD-HEP-99-6 (1999), tiep-ph/9903294. 



[15] D. Atwood, S. Bar-Shalom and A. Soni, Brookhaven preprint BNL-HET-99/11, hep- 



ph/9903538 



[16] S. Cullen and M. Perelstein, SLAC preprint SLAC-PUB-8084 (1999), hep-ph/03422. 



[17] T.G. Rizzo, SLAC preprint SLAC-PUB-8036 (1998), hep-ph/9901209 



[18] J.L. Hewett, SLAC preprint SLAC-PUB-8001 (1998), |hcp-ph/98H356| . 

[19] P. Mathews, S. Raychaudhuri and K. Sridhar, TIFR preprint TIFR-TH/99-13 (1999), 
frep-ph/9904232 . 



[20] P. Mathews, S. Raychaudhuri and K. Sridhar, TIFR preprint TIFR-TH/98-48 (1998), 
tiep-ph/9811501| (to appear in Phys. Lett. B). 

[21] P. Mathews, S. Raychaudhuri and K. Sridhar, TIFR preprint TIFR-TH/98-51 (1998), 
ficp-ph/9812486| , (to appear in Phys. Lett. B); T.G. Rizzo, SLAC preprint SLAC-PUB- 



7986 (1999), |hep-ph/9902273 



[22] K. Agashe and N.G. Deshpande, Oregon U. preprint OITS-669 (1999), frep-ph/ 9902263 



[23] C. Balasz, et al, Michigan State U. preprint MSUHEP-90105 (1999), |hcp-ph/9904220| . 



[24] T.G. Rizzo, SLAC preprint SLAC-PUB-8071 (1999), |hciTph/9903475 
[25] H.L. Lai, et al., Phys. Rev. D55, 1280 (1997) . 

[26] R. Hamberg, W.L. van Neerven and T. Matsuura, Nucl. Phys. B359, 343 (1991) 



[27] B. Abbott, et al, D0 Collaboration, Fermilab preprint FERMILAB-PUB-98/391-E, hcp- 



ex/9812010 



[28] F. Abe, et al, CDF Collaboration, Phys. Rev. D59, 052001 (1999) 



14 



